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ABSTRACT 

Near-infrared polarimetric imaging observations toward the Galactic center 
have been carried out to examine the efficiency and wavelength dependence of 
interstellar polarization. A total area of about 5.7 deg 2 is covered in the J, H, 
and Kg bands. We examined the polarization efficiency, defined as the ratio 
of degree of polarization to color excess. The interstellar medium between the 
Galactic center and us shows the polarization efficiency lower than that in the 
Galactic disk by a factor of three. Moreover we investigated the spatial variation 
of the polarization efficiency by comparing it with those of color excess, degree 
of polarization, and position angle. The spatial variations of color excess and 
degree of polarization depend on the Galactic latitude, while the polarization 
efficiency varies independently of the Galactic structure. Position angles are 
nearly parallel to the Galactic plane, indicating the longitudinal magnetic field 
configuration between the Galactic center and us. The polarization efficiency 
anticorrelates with dispersions of position angles. The low polarization efficiency 
and its spatial variation can be explained by the differences of the magnetic 
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field directions along the line-of-sight. From the lower polarization efficiency, we 
suggest a higher strength of a random component relative to a uniform component 
of the magnetic field between the Galactic center and us. We also derived the 
ratios of degree of polarization Ph/pj = 0.581 ± 0.004 and Pk s /ph = 0.620 ± 
0.002. The power law indices of the wavelength dependence of polarization are 
(3jh = 2.08 ± 0.02 and (3hk s = 1-76 ± 0.01. Therefore the wavelength dependence 
of interstellar polarization exhibits flattening toward longer wavelengths in the 
range of 1.25—2.14 /jm. The flattening would be caused by aligned large-size 
dust grains. 

Subject headings: polarization — dust, extinction — ISM: magnetic fields - 
Galaxy: center — infrared: stars 



1. INTRODUCTION 



Observations of interstellar linear polarization (hereafter just interstellar polarization) 
give information about the magnetic field and properties of polarizing dust grains. Interstel- 
lar polarization is caused by n on-spherical dust gra ins aligned by a magnetic field (dichroic 
extinction; see e.g., reviews by lLazarianll2003l 120071 ). 



Pos ition angles of polar i zation yield the directions of magnetic fields on the plane-of- 
the-sky. iMathewson fc Fordl (119701 ) compiled the polarization data of nearly 7000 stars and 
showed the distribution of position angles in the Galactic coordinates. The distribution of 
position angles for stars beyond 1 kpc of the Sun traces the structure of the large-scale 
Galactic magnetic field which runs almost parallel to the spiral arms, while that for stars 
within 600 pc traces the structure of the local magnetic field. The structure of magnetic 
field is less simple at / ~ 40°, 80°, 260°, and away from the Galactic plane such as the active 
star forming regions Taurus, Perseus, Ophiuchus, and Orion, where position angles deviate 
from the large-scale l ongitudinal p attern. The local magnetic field points toward / ~ 80° and 
away from I ~ 260° (jHeileslll976l ). From the direction of I ~ 40°, a loop structure extends 
toward the north Galactic pole. Toward the loop structure and the star forming regions, the 
magnetic field is obviously perturbed from the uniform, large-scale structure. These indicate 
the existence of a random component of the magnetic field on a small-scale in addition to a 
uniform component on a large-scale. 

The random component of the magnetic field can cause a dispersion in the measured 
position angles and a decrease of the polarization efficiency. The polarization efficiency is 
defined as a ratio of polarization degree to extinction such as P\/t x and p\ 2 /E(\i — A 2 ). If 
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magnetic fields are tangled along the line-of-sight, degree of polarization p\ does not build 
up as much as it would in a uniform magnetic field. The extinction (as measured by t\ or 
E(\i — A 2 )) along a line-of-sight is determined only by the total column of dust, so that 
the polarization efficiency decreases (depolarization; iMartinl 11974 ) . Thus, the polarization 
efficiency is a useful measure to probe the random component of the magnetic field. 



Serkowski et al.l (119751 ) analyzed a sample of 180 nearby stars that were observed po- 



larimetrically in the U, B, V, and R bands. They investigated the relation between color 



excess E(B — V) and the maximum polarization p E 
upper limit for the polarization efficiency, 



at the wavelength A max , and found an 



Pu 



JE(B -V)< 9.0%/mag. 

The polarization efficiency changes from line-of-sight to line-of-sight below the upper limit. 

Observations at optical wavelengths are limited to regions with small extinction (typi- 
cally Ay < 5 mag). The be havior of the p olari zation efficiency at large extinction (up to Ay 
~ 100 mag) was studied by IJoned (Il989l ) and IJones et al.l (119921 ) using near-infrared (NIR) 
wavelengths. They compiled K band (2.2 //m) polarimetric measurements for about 100 
sources at various locations and extinctions. The relation between polarization and extinc- 
tion was modeled by assuming that interstellar polarization depends only on the geometry 
of the uniform and random components of the magnetic field. The model fits the observed 
relation between polarization and extinction well when equipartition of the energy density 
holds between the uniform and random components of the magnetic field. 

Studies on the polarization efficiency of the diffuse interstellar mediu m (ISM) hay e 
been carried out using the comprehensive compilation of polarization data of iHeilea (120001 ) . 
Based on the polarization data for abo ut 5,500 stars distrib uted over the entire sky and 
mostly located at distances of d < 4 kpc, iFosalba et al.l (120021 ) found a nearly linear growth 
of average polarization degree with extinction up to E (B — V) ~ 1 mag, but noted that the 
polarization efficiency is much lower than what is expected from completely aligned grains 
in a uniform magne tic field. They found th e polarization efficiency that was 1/3 that of the 
maximum found by lSerkowski et al.l (Il975l ). They explained it by depolarization due to the 
random component of the magnetic field, and estimated the magnetic field strength ratio of 
the uniform to the random component, B u /B r , to be about 0.8, where B u and B r are the 
strengths of the uniform and random components. 

In the innermost region of the Gal a xy, only a few pola ri metric observations haye been 
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done (e.g., 

measured polarization for a few dozen of highly reddened (H — K < 3 mag) stars in the area 
at / ~ 0°, 20°, and 30° in the Galactic plane, and noted that the polarization efficiency is 
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lower by a factor of about four than t hat in the sol a r neig hborhood. Based on polarization 
measurements for 127 reddened stars, I Creese et al.l (119951 ) reported that greater extinction 
results in increased polarization, but the increase is smaller than expected. They concluded 
that the polarization efficiency is lower than that in the solar neighborhood. Their samples 
are not highly extincted (H — K < 0.7 mag), because they were selected from an I band 
objective prism survey and were therefore relatively sparse and shallow. 

The properties of polarizing dust grains can be examined by the wavelength depen- 
dence of interstellar polarization. Polarization p\ shows a convex curvature with a peak 
jPmax, typically occurri ng at around A max = 0.55 /xm, and a wing toward NIR wavelengths. 
ISerkowski et al.l ( 19751 ) made a determination of the wavelength dependence of polarization, 
empirically establishing "Serkowski's law", 



Px/Pmsx = exp[-1.151n 2 (A n 
The slope of the wing is represented by a power law as 

Px oc \~ f3 



<A)]. 



1992 


)■ 


Creese et al. 


(1995) 



for about 10 reddened stars with respect to a power law. However, more samples should be 
needed to confirm this suggestion. 

Now polarimetric imaging observations using several hundred thousand stars lying in 
the Galactic bulge as background sources enable us to measure polarization degrees, position 
angles, and color excess finely across the area and deeply to the line-of-sight. Therefore, we 
carried out NIR polarimetric imaging observations for a wide field (~ 3° x 2°) toward the 
Galactic center (GC). From the polarization efficiency and dispersions of position angles, 
we discuss the uniform and random components of the magnetic field along the line-of- 
sight on the Galactic scale. Furthermore, we verify whether the wavelength dependence of 
polarization at NIR is really represented by a power law or not, and examine what types 
of dust grains cause the NIR polarization. In section 2, we describe observations and data 
reduction. In section 3, we present results for color excess, polarization, the polarization 
efficiency, and the wavelength dependence of polarization in our sample. In section 4, we 
discuss the results as they relate to the magnetic field and polarizing grain properties toward 
the GC. 
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2. OBSERVATIONS AND DATA REDUCTION 

Polarimetric imaging observations toward the GC have been carried out with a NIR 
polarimetric camera installed on the IRSF (InfraRed Survey Facility) 1.4 m telescope at 
the South African Astronomical Observatory (SAAO) in Sutherland. The NIR polarimet- 



Kandori et al. 


2006 


and NIR camera SIR- 


for Unbiased Survey: 


Nagashima et al. 


1999; 



Nagayama et al.ll2003l ). The camera is equipped with three 1024 pixel x 1024 pixel HAWAII 
arrays. This enables simultaneous observations in the J (central wavelength Xj = 1.25 jum), 
H (Xh = 1-63 jum), and K s (A# s = 2.14 //m) bands by splitting the beam into the three 
wavelengths with two dichroic mirrors. The image scale of the arrays is 0"45 pixel -1 , yielding 
a field of view of 460" x 460". 

From 2006 to 2009, we have observed 459 fields toward the GC, and the total area 
covered is about 5.7 deg 2 (see Fig. [Q. The centers of fields were set at intervals of 400". 
We obtained 10 dithered frames on the circle with a radius of 20", yielding an effective field 
of view of about 420" and overlaps between adjacent fields with a size of about 420" x 20". 
We performed 10 s exposures at four wave plate angles (0°, 45°, 22?5, and 67?5), resulting 
in a total exposure time of 100 s per wave plate angle for each field. Our observations were 
carried out under stable sky conditions on photometric nights. The seeing was typically 1"3, 
1"2, and 1"1 (FWHM) in the J, H, and Ks bands, respectively. To make median sky frames 
we observed one of two sparse stellar fields (I = -2?933, b = 7?100; I = 4?525, b = -12?838) 
for each one or two field(s), that is, as frequently as every about 10 or 20 minutes. Twilight 
flat frames were obtained before and after the observations. Dark frames were obtained at 



the end of the nights. The polarimetric standard star R CrA No. 88 (IWhittet et al.lll992l ) 
was observed 15 times through the observing runs, with 1.6 or 2 s exposures at each wave 
plate angle at 10 dithered positions. 

We applied the standard procedures of NIR array image reduction, including dark- 
current subtraction, flat-fielding, sky subtraction, and frame combination using the IRAF 
(Image Reduction and Analysis Facility)© software package. After subtraction of an averaged 
dark frame, each frame was divided by a normalized flat frame. Then the thermal emission 
pattern, the fringe pattern due to OH emission, and the reset-anomaly slope pattern of the 
HAWAII arrays were subtracted from each frame with a median sky frame. This subtraction 



6 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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cannot be done adequately in the case that the intensity of OH emission or temperature 
abruptly change, so that we observed the field again. Finally we obtained images by com- 
bining 10 frames at each wave plate angle, and Stokes I images by combining 10 x 4 = 40 
frames. 

Photometry of point sources was performed using the DAOPHOT package in IRAF. 
We used the DAOFIND task to detect point sources on Stokes I images. Since the observed 
fields are highly crowded stellar fields, we obtained positions and magnitudes of detected 
sources on Stokes I images from PSF-fitting photometry using the ALLSTAR task. A 
model PSF was constructed from bright and isolated sources whose numbers were typically 
20—60, 50—90, and 60—90 on each image in the J, H, and K$ bands, respectively. Positions 
on ce lestial coordinate syst ems were calculated referring to the 2MASS Point Source Cata- 
log0 ( Skrutskie et al. 2006 ). The positional accuracy was estimated to be about 0'.'03. For 



photometric calibration, we compared magnitudes for the detected sources with those for 
the 2MASS Point Source Catalog sources in the observed field. In this comparison magni- 
tude transformations from the 2MASS system to the IRSF system were applied (Y. Nakajima 
2007, private communication). The photometric accuracy is about 0.03 mag in all the JHKs 
bands. 

In order to estimate the effect of confusion on the PSF-fitting photometry in crowded 
fields, we computed the completeness of the PSF-fitting photometry. We added 900 artificial 
sources with the brightness corresponding to the limiting magnitudes, J ~ 14.0 mag, H ~ 
13.4 mag, and K$ ~ 12.5 mag (see the last paragraph in this section), on each Stokes I image 
in a reticular pattern at intervals of 14" (30 x 30 sources) using the ADDSTAR task, and 
performed source detection and PSF-fitting photometry in the same manner as above. The 
recovery rates of the added artificial sources and median differences of magnitudes between 
the recovered and added artificial sources were calculated in each field. The recovery rates 
were 98.5%, 95.9%, and 95.6% on average, and the means of the absolute median differences 
of magnitudes were 0.005 mag, 0.011 mag, and 0.009 mag in the J, H, and Ks bands, 
respectively. 

In order to calculate the Stokes parameters for sources, we performed detection and 
aperture photometry of sources at each wave plate angle using the DAOFIND and PHOT 
tasks, and obtained flux and positions. We compared aperture photometry with PSF-fitting 
photometry by calculating normalized Stokes parameters Q/I and U/I for duplicate sources 
in overlapping regions of adjacent fields observed under different observing conditions. As a 



7 We related pixel coordinates to celestial coordinates using the OPM software, wh ich is compile d by Dr. 
N. Matsunaga and based on the Optimistic Pattern Matching algorithm proposed by iTaburl ( 2007 ). 
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result, aperture photometry gives a better result than PSF-fitting photometry which shows 
systematic offsets. The center of aperture was determined as the centroid of each source 
independently on images at each wave plate angle. We adopted the aperture size of 1.5 x 
FWHM measured on Stokes I images. The aperture size was chosen after search of 1.0, 1.5, 
and 2.0 x FWHM. The normalized Stokes parameters were best matched for a factor of 1.5 
with each other between overlapping regions (see below for the comparison of normalized 
Stokes parameters in overlapping regions). 

We also computed the completeness of the aperture photometry in a similar manner 
to the PSF-fitting photometry. We added 900 artificial sources with the brightness corre- 
sponding to the limiting magnitudes on each image at each wave plate angle in a reticular 
pattern at intervals of 14" (30 x 30 sources), and performed source detection and aperture 
photometry as described above. We calculated the recovery rates of the added artificial 
sources and standard deviations of median flux of the recovered artificial sources among all 
wave plate angles in each field. The average recovery rates were 97.5%, 94.4%, and 94.3%, 
and the standard deviations of the median flux were 0.2%, 0.3%, and 0.2% on average in the 
J, H, and K s bands, respectively. 

We merged sources on Stokes / image from PSF-fitting photometry and sources on 
images at each wave plate angle from aperture photometry. Positions of sources obtained 
with PSF-fitting photometry on Stokes I images were used as the reference positions for 
merging. We matched sources on Stokes / image and images at each wave plate angle to the 
reference positions within a 1 pixel radius. In crowded fields, PSF-fitting photometry yields 
more precise positions of sources than DAOFIND. We mitigated source confusion by using 
the precise reference positions for matching. Based on flux of sources at each wave plate 
angle, we calculated Stokes parameters /, Q, U, and their statistical errors (calculated from 
noise of signal, sky background, dark current, and readout). The Stokes parameter I was 
calculated from the flux on images at each wave plate angle, as well as the Stokes parameters 
Q and U, not from the flux from the PSF-fitting photometry on Stokes / images, in order to 
minimize error due to variation of PSF among images at each wave plate angle and Stokes 
I images. These Stokes parameters also have systematic errors originating from change of 
atmospheric conditions (seeing and transparency of atmosphere) at each wave plate angle 
for each field. To estimate these systematic errors we have made a comparison of normalized 
Stokes parameters Q/I and U/I of the same sources in overlapping regions between adjacent 
fields with sizes of about 420" x 20". We defined the same sources as the sources having 
closest positions within 1". Using the sources whose statistical errors of Q/I and U/I were 
less than 1%, average differences of Q/I and U/I were calculated in each overlapping region. 
The numbers of the same sources in each overlapping region were typically 10—50, 40—100, 
and 40—100 in the J, H, and Ks bands, respectively. The means of the average differences 
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of Q/I and U/I in overlapping regions were 0.33% and 0.35% in the J band, 0.26% and 
0.27% in the H band, and 0.25% and 0.26% in the K s band. We adopted these values 
for the systematic errors of the normalized Stokes parameters Q/I and U/I, and computed 
the total errors by combining the estimated systematic errors with the statistical errors in 
quadrature. Then the observed degree of polarization p Q b s , position angle 9 and their errors 
Sp and 89 were derived. 9 is defined as the angle between ^-vector of polarization and the 
direction of north celestial pole and increasing to the east (—90° < 9 < 90°). To correct 
noise biasing we calculated degree of polarization p using the following formula: 

P = VPobs 2 - Sp 2 



( iWardle fc KronbergJ 11974 ; Oarke fc Steward Il986l ). We regarded sources with p ^ s < Sp as 
unpolarized (p = 0) sources. All the data were calibrated for the polarization efficiency of 
the wave plate and polarizer (95.5%, 96.3%, and 98.5% in the J, H, and Kg bands; see 
Kandori et al1l2006h . 



We checked our analysis b y comparing our pola rimetry of the polarimetric standard 
star R CrA No. 88 with that by lWhittet et al.l (119921 ). The values p and 9 were derived by 
averaging normalized Stokes parameters Q/I and U/I obtained from 15 observations. Their 
errors 5p and 69 were determined from the standard deviations of the means of Q/I and 
U/I. As shown in T able (U we confirmed that our polarimetry is consistent with that of 
Whittet et al.l (119921 ) within the errors in all the bands. 



For an additional check of our analysis, we made a comparison of p and 9 of the same 
sources in overlapping regions in a similar way to the comparison of Q/I and U/I. Using 
the same sources with Sp < 1% and 59 < 10°, we calculated differences of p and 9. The 
numbers of the same sources in each overlapping region were typically 10—30, 20—80, and 
10—70 in the J, H, and Ks bands, respectively. Figure [2] shows differences of p and 9 as 
a function of means of p for the sources. In all the range of the measured polarization, the 
absolute differences of p are mostly less than 1%, and those of 9 are mostly less than 10°; 
the standard deviations of the differences of p are 0.79%, 0.55%, and 0.49%, and those of 9 
are 6?9, 5?3, and 6?3 in the J, H, and Ks bands, respectively 

Source detections at J, H, and K s were merged into a single source record using the 
positions in each band. First the Ks sources were taken as seed detections; then the Ks-H 
pairwise matching was done, and this was followed by the Ks-J pairwise matching. Then 
the H sources which were not matched in the previous process were taken as seeds; then 
the H-J pairwise matching was done. The match was acceptable if the source separation 
was less than 1". We adopted the coordinates for the longest wavelength in the matching 
as the source coordinates. There are a total of 3,539,087 sources: 1,536,017, 2,979,994, and 
3,190,511 sources detected in the J, H, and Kg bands, respectively. The limiting magnitudes 
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of our survey, defined as the level at which 5p < 1%, are J ~ 14.0 mag, H ~ 13.4 mag, and 
K s ~ 12.5 mag. There are 234,121 sources with 5pj < 1%, 541,990 sources with 5pu < 1%, 
and 558,647 sources with 6px s < 1%- 



3. RESULTS 

3.1. MK CLASSIFICATIONS AND LOCATIONS 

A color-color diagram for the sources is shown in Figure |3j The figure includes 165,858 
sources that have 5p < 1% in all the bands. There are three distinct populations; the majority 
of the sources fall in a feature extending parallel to the reddening vector from the locus of 
giants, while two weak concentrations of the sources are seen around (H — K s , J — H)~(0.1, 
0.3) mag and ~ ( 0-2, .7) mag along the loci of dwarfs and giants. Based on the model by 



Wainscoat et al.l (119921 ). we estimate what kind of sources can be detected in all the bands 
within the limiting magnitudes (J ~ 14.0 mag, H ~ 13.4 mag, and Ks ~ 12.5 mag). The 
model predicts the numbers of sources expected for each MK classification (spectral type and 
luminosity class) at each distance from the Sun. Toward the GC, we would expect that most 
of the detectable sources are K/M giants located in the Galactic bulge. They are heavily 
reddened due to the large amount of the intervening dust, corresponding to the extended 
feature in Figure [3j Among the sources located in the Galactic disk (< 4 kpc from the Sun), 
A/F dwarfs and G/K giants are mainly detected. The two concentrations around the loci of 
dwarfs and giants predominantly consist of these dwarfs and giants with small extinction. 

Out of 234,121 sources with 5pj < 1%, 541,990 sources with 5p H < 1%, and 558,647 
sources with 5pK s < 1%, those having H — Ks colors amount to 196,651, 512,030, and 
544,675 sources in the J, H, and K s bands, respectively. We divide these sources into two 
groups at H — Ks ~ 0.4 mag, where a saddle of the extended feature of giants exists; the 
disk sources have H — Ks < 0.4 mag and bulge sources have H — Ks > 0.4 mag. The 
former are mostly A/F dwarfs and G/K giant located in the Galactic disk, while the latter 
are mostly K/M giants located in the Galactic bulge. The numbers of the disk and bulge 
sources amount to 58,007 and 138,644 with 5pj < 1%, 54,102 and 457,928 with 5p H < 1%, 
and 31,448 and 513,227 with 5px s < 1%. Hereafter we call these sources as the disk and 
bulge sources. 
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3.2. COLOR EXCESS 



We calculate color excess for the disk and bulge sources using the equation: 

E(H - K s ) = {H- K s ) - ((H - K s ) ). 

The mean of the intrinsic colors of sources, ((H — K$)o), is computed as follows. First we 
compute the numbers of sources e xpected in each b a nd for each MK classification in the disk 
and bulge based on the model by IWainscoat et al.l (119921 ) under the following criteria. The 
criteria for the disk and bulge sources are detection within the limiting magnitudes in each 
band, that is, detection with J < 14.0 mag, H < 13.4 mag, and Ks < 12.5 mag. Additional 
criterion for the disk sources is H — K$ < 0.4 mag, and that for the bulge sources is H — Ks > 
0.4 mag. We compute {(H — Ks)o) by averaging the intrinsic colors of sources, (H — Ks)co, 
using the equation: 

^{(F-if s ) co xN c } 
((H - K s )o) = ^ ^ , 



where the sum is over the 29 spectral classes in the IWainscoat et al.l ( 19921 ) model with 
luminosity classes of IH and V (see their Table 2), Nc is the total number of sources of class 
C predicted by th e mo del, and the intr i nsic c olors, (H — Ks)co, for each class are taken from 
Koornneeil (119831 ) and lBessell fc Brettl (119881 ) . For the uncertainty of ((H — K s )o), we take 
the standard deviation of the intrinsic colors of sources from an equation: 



a((H-K s )o) 



- A»co - ((H - Ksk)} 2 x No] 



((H — Ks)o) and <y((H — Ks)o) in an area with a size of 10' x 10' are computed in each 
direction of the disk and bulge sources. From line-of-sight to line-of-sight, ((H — Ks)o) 
and <t{{H — Ks)o) vary, because of the spatial variations of number density of sources with 
each MK classification in the Wainscoat's model. As for the disk sources, ((H — Ks)o) and 
<r((if — Ks)o) are computed to be 0.06—0.08 mag and 0.05 mag in the J band, 0.07—0.08 
mag and 0.04—0.05 mag in the H band, and 0.08—0.10 mag and 0.04—0.05 mag in the Ks 
band for lines of sight toward each disk source. As for the bulge sources, ((H — K s )o) and 
<r((if — K s )o) are computed to be 0.13—0.19 mag and 0.05 mag in the J band, 0.16—0.21 
mag and 0.04—0.08 mag in the H band, and 0.17—0.22 mag and 0.04—0.06 mag in the Ks 
band for lines of sight toward each bulge source. 
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The color excess E(H — Kg) for a disk source is not so large compared with its error, 
and therefore cannot be accurately determined. We adopt only means of E(H — Kg) for 
the disk sources for further discussions. From the means of H — Kg (0.20 mag, 0.21 mag, 
and 0.22 mag) and average intrinsic color of the disk sources (0.07 mag, 0.08 mag, and 0.09 
mag), the means of E(H — Kg) are calculated to be 0.13 mag, 0.13 mag, and 0.13 mag for 
the disk sources in the J, H, and Kg bands, respectively (dash-dotted lines in Fig. H]). 

Contrary to E(H — Kg) for a disk source, E(H — Kg) for a bulge source is so large that 
it can be accurately determined. We show histograms of E(H — Kg) for the bulge sources 
in Figure HJ The distributions of E(H — Kg) peak at about 0.5 mag and have tails toward 
large color excess in all the bands. The bulge sources in the H and Kg bands trace larger 
color excess than those in the J band because of the wavelength dependence of interstellar 
extinction. 

We present a map of E(H — Kg) for the bulge sources in the Kg band as follows. The 
observed area is divided into cells with a size of 2' x 2'. Means of E(H — Kg) for the bulge 
sources in the Kg band are calculated in each cell. The size of cells is determined to contain 
as many sources in each cell as possible without degradation of the angular resolution. The 
resultant map of the mean (E(H — Kg)) is shown in Figure |5j The map shows the spatial 
variation of interstellar extinction depending on the Galactic latitude. The cells with large 
values (E(H — Kg)) > 2 mag mainly concentrate to the Galactic plane (|6| < 0?5), showing 
clumpy and filamentary structures. 

3.3. DEGREE OF POLARIZATION AND POSITION ANGLE 

We plot a Kg band polarization vector map for the disk and bulge sources in Figure [6j 
Zooming the electronic edition of the map makes the .E-vectors of polarization for sources 
legible. Most of the sources show the 22- vectors of polarization nearly parallel to the Galactic 
plane, while there exist some deviations. These deviations of .E-vectors are seen in the 
regions where a relatively few sources with large color excess are detected. Most of the 
detected sources with large color excess could be in/behind the nearby dense clouds. The 
deviated .E-vectors would reflect the local magnetic field directions in nearby dense clouds. 

Figure [7] shows the relations between position angles and colors for the disk and bulge 
sources with 59 < 10° in each band. There are two distinct populations in blue (H — Kg < 
0.4 mag) and red (H — Kg > 0.4 mag) colors, corresponding to the disk and bulge sources, 
respectively (see §3.11 and Fig. [3]). The disk sources are prominent in the J band, but less 
prominent in the H and Kg bands. The numbers of the disk and bulge sources amount 
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to 38,425 and 118,045 with 66 j < 10°, 20,858 and 355,124 with 86 H < 10°, and 6,086 and 
323,603 with 59 k ? < 10° ■ The mea ns of 9 clearly differ between the two populations as 
also noted by lKobayashi et al.l ( 119831 ) and iNishiyama et al.l (j2009al ); the means of 9 are 9?3, 
10?3, and 11? 4 for the bulge sources, while — 1?6, — 2?1, and 0?0 for the disk sources in the 
J, H, and K$ bands, respectively. Most of the bulge sources have 9 nearly parallel to the 
Galactic plane (~27°), but slightly rotated westward. Position angles rotate more for the 
disk sources. 

In a similar manner as E(H — K$), we present a map of 9k s for the bulge sources with 
S0k s < 10 ° in Figure HI Most of cells show the values 0° < (9 Ks ) < 20°, which indicate that 
the magnetic fields between the GC and us are longitudinal on average. These values are 
close t o the average position a ngle which traces the magnetic field configuration in the GC 



(16.0°; INishiyama et al.ll2009a ). suggesting that the Galactic magnetic field running nearly 



Novak et al. 


2003; 


Chuss et al. 


2003) 



In the relations between polarization degrees and colors for the disk and bulge sources 
(Fig. [H]), two distinct populations corresponding to the disk and bulge sources can be also 
seen. The disk sources have average polarization degrees of 2.1%, 1.2%, and 0.8% in the J, 
H, and Ks bands, respectively. Meanwhile, the bulge sources extend red-ward and show 
correlations between p and H — Ks in all the bands; p increases with increasing H — Ks- The 
slopes of the correlations correspond to the polarization efficiency ( §3.4p . and the difference of 
the slopes between the bands is due to the wavelength dependence of interstellar polarization 
(ffl). 

In Figure [T0| we show a map of the mean (pk s )- In a similar way to the spatial variation 
of E(H — Ks) (Fig. [5]), that of px s is dependent on the Galactic latitude. Most of cells 
in which (pk s ) exceeds 5% are close to the Galactic plane, while in cells at higher Galactic 
latitude (pk s ) is only 1—2% or less. However, the spatial variations of E(H — Ks) and p 
do not completely coincide with each other. This is more obvious in the form of the spatial 
variation of the polarization efficiency ( §3.41) . 



3.4. THE POLARIZATION EFFICIENCY 

Starlight suffers both extinction and polarization in the passage through the intervening 
ISM. Starlight from the disk sources passes through the ISM in the disk (< 4 kpc from the 
Sun), whereas starlight from the bulge sources passes through the ISM between the GC and 
us (i.e, the ISM in the disk and bulge; see §3.1j) . From p/E(H — Ks) for the disk and bulge 
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sources, we examine the polarization efficiency of the ISM in the disk and that between the 
GC and us. 

Combining means of E(H — Kg) and p, we calculate the means (p)/(E(H — Kg)) for 
the disk sources to be 16.2% / mag, 9.2% / mag, and 6.2% / mag in the J, H, and K s 
bands, respectively (dash-dotted lines in Fig. [TT1) . We show histograms of p/E(H — Kg) for 
the bulge sources in Figure [TTJ The means of p/E(H — Kg) for the bulge sources are 6.3% 
/ mag, 3.7% / mag, and 2.4% / mag in the J, H, and Kg bands, respectively (shown as 
arrows in the histograms). The average polarization efficiency of the ISM between the GC 
and us is considerably lower than that of the ISM in the disk, by a factor of about three. 

We make a comparison between observed polarization efficiency and estimated upper 
limits [pj/E(H~K s ) = 25.0% / mag, p H /E(H-K s ) = 14.5% / mag, and p Ks /E(H-K s ) = 
9.0% / mag], which are e stimated by extending the upper limit \p max /E(B— V) = 9.0% / mag] 
at optical wavelengths (ISerkowski et al.lll975l ) to NIR wavelengths as follows. We convert 
color e xcess E(B — V) to E(H — Kg) by assuming that E(H — K) of the interstellar extinction 
law bv lRieke fc Lebofskvi fjl985h is identical to E(H - Kg); E(B - V)/E(H - Kg) = 5.14. 
As for degree of polarization, first, Pp ax is converted to pk (th e value at 2.2 iim) using the 
average ratio(p max /pif) = 5.4 (|Joneall989l ; IWilking et al.lll980l ). Then px s are extrapolated 
from pk following a power law p\ oc A" 1 ' 76 , and pj and pu are calculated from px s using 
Ph/pj = 0.581 and Pk s /ph = 0.620 ( §3.5j) . The average polarization efficiency of the ISM 
in the disk is about two-thirds of the upper limits. Moreover, that of the ISM between the 



GC an d us is no more than about a quarter of the estimated upper limits. iKobayashi et al. 



(119831 ) also suggested that the polarization efficiency toward the GC is considerably lower 
than that obtained in the solar neighborhood based on K band polarimetry toward the GC 
(20' x 20'). 

The standard deviations of polarization efficiency a(p/E(H — Kg)), which are larger 
than the average errors of p/E(H — Kg), show that the polarization efficiency has spatial 
variation. As for E(H — Kg), 9k s , and px s , we present a map of px s /E(H — Kg) for the 
bulge sources in Figure [121 It shows large variation from line-of-sight to line-of-sight in a 
range of about 1 to 5% / mag, but does not depend on the Galactic longitude and latitude, 
nor on the Galactic structure. 

Dispersions of position angles cf{9k s ) anticorrelate with {pK s /E(H — Kg)). We calculate 
a{9 Ks ) for the bulge sources with 59 k s < 10° in each cell and show its map in Figure [131 
Some cells show <t(6k s ) significantly larger than the average error of 0k s for the bulge sources 
(5?7). In comparison between the spatial variations of pk s /E(H — Kg) and <t(9k s ) (Figs. [T21 
and [T3|) . we can see a tendency that the larger <t(9k s ), the lower (pk s /E(H — Kg)) , and vice 
versa. This tendency is shown in the relation between (pk s /E(H — Kg)) and <t(9k s ) (Fig. 
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IT4"1) . The medians of (pk s /E(H — Kg)) in each bin of ct(6k s ) (2° in width) show the highest 
value at the bin of a{6 Ks ) = 2—4°, decrease with increasing a(9 Ks ), and then become almost 
flat toward bins of larger <j(9k s )- 



3.5. THE WAVELENGTH DEPENDENCE OF POLARIZATION 

We here examine the wavelength dependence of polarization. Out of the detected 
sources, we select 3,651 sources that have 5p < 1% and p > 10 5p in all the bands. Here- 
after we use these sources in this section. 3,647 of these sources have H — K$ > 0.4 mag, 
and therefore they are bulge sources ( §3.ip . Figure [TBI shows the spatial distribution of the 
sources. The distribution of these bulge sources is little clumpy on regions where extinction 
is moderate (not especially small and large) but the sources are seen in diverse area in /, b 
(see also Fig. [5j. In Figure [TBI the correlations pj vs. pu and pu vs. px s are good with linear 
regressions of (p H /pj) = 0.581 ± 0.004 and (pk s /ph) = 0.620 ± 0.002 (Table[2]). The scatter 
around the best-fitting lines are 0.076 and 0.047 for Ph/pj and Pk s /Ph, respectively. The 
scatter can be explained by the observational error in p\\ the average errors are (S(ph/pj)) 
= 0.055 and (6(pk s /ph)) = 0.063. Thus the wavelength dependence of polarization does not 
change significantly from line-of-sight to line-of-sight in our sample. 

Assuming a power law (p x oc \~@), and using the equations: 

\n(p H /pj) 



0. 



JH 



3 



HK S 



Mpk s /ph) 



\n(X Ks /X H ) 

we calculate indices (3jh and (5hk s for the sources. Histograms of (3jh and (5hk s are shown 
in Figure [T71 The means {(3jh) and (Phk s ) for the sources are 2.08 ± 0.02 and 1.76 ± 0.01, 
respectively (Table [2]). The errors of the means and the average errors of (3jh and (3hk s for 
the sources are derived in a similar manner to Ph/pj and Vk s Ivh- Although these values 



are not inconsistent with the empirical values of 1.6—2.0 (lNagatalll990l : Martin &: Whittet 



19901 ; iMartin et al.lll992l ). {(3jh) is larger than (/3hk s )- The degree of polarization decreases 
more slowly than a power law as the wavelength becomes longer from 1.25 to 2.14 /im; the 
wavelength dependence of polarization appears to flatten toward longer wavelengths. 
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4. DISCUSSION 

4.1. LOW POLARIZATION EFFICIENCY AND ITS SPATIAL VARIATION 



As for the polarization efficiency of the ISM between the GC and us, we revealed (a) 
low efficiency compared to that of the ISM in the disk (Fig. [TT1) . (b) the spatial variation 
throughout the observed area (Fig. [T2|) . and (c) the anticorrelation with the dispersions 
of position angles (Fig. To explain our results, we discuss the polarization efficiency 

in relation to two factors: (1) the polarizing grain properties and (2) the magnetic field 
direction. 

The polarizing grain properties such as their shape, size distribution, and composition 
could affect the polarization efficiency. However, the wavelength dependence of polarization 
shows little or no spatial variation in our sample ( §3.5p . This suggests that the polarizing 
grain properties are almost uniform, and that the result of (b) cannot be explained by the 
factor (1), the polarizing grain properties. 

Superposition of the ISM with different magnetic field directions along the line-of-sight 
could affect the polarization efficiency. The lines-of-sight towards the disk and bulge sources 
crosses multiple ISM with a range of physical conditions. The interstellar magnetic field 



consists of the uniform and random components (jHeileal 19871 . Il996f ) . The uniform component 



( Heiles 


1996; 


Han 


2009) 



as a part of ISM; in each segment, both the magnetic field direction and the degree of grain 
alignment are constant) can deviate from the direction of the Galactic magnetic field due 
to the presence of a random local component. Due to the differences of the magnetic field 
directions, superposition of the ISM along the line-of-sight lowers the polarization efficiency 
(depolarization). Since the ISM between the GC and us generally consists of more ISM 
with different magnetic field directions, depolarization would be larger and the polarization 
efficiency should be lower than the ISM in the disk (result of (a)). To explain the observed 
dispersions of position angles (Fig. [13]), nonuniform str uctures of the magnetic fie ld and/or 



density with a size of less than 2' (cell size) are needed. iGosling et al.l (120061 . 120091 ) detected 
such nonuniform density distribution with a size of 5—15". Larger (smaller) differences of 
the magnetic field directions would cause lower (higher) polarization efficiency and a larger 
(smaller) dispersion of position angles (results of (b) and (c)). Of these two factors, only the 
factor (2), the differences of the magnetic field directions along the line-of-sight, can explain 
our results. 
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4.2. THE MAGNETIC FIELD STRENGTH RATIO OF THE RANDOM TO 

THE UNIFORM COMPONENT 



We discuss the magnetic field strength ratio of the random to the uniform component 
based on the observed re lation between extinction and degree of polarization. For the re- 
lation, I Jones et al.l (119921 ) constructed two models depending on the geometry of magnetic 
fields along the line-of-sight. One is the two-component model, and the other is the wave 
model. 

In the two-component model, the magnetic field direction is determined by a combina- 
tion of the uniform and random components in each optical depth length (segment); in each 
length Atk s = 0.1, the random component of the magnetic field decorrelates. A segment 
corresponds to a part of the diffuse ISM (with a typical len gth of a few tens of pc) or a 
dense cloud (a fraction of a pc). Fitting the model to the data, IJones et al.l (119921 ) concluded 
that the uniform and random components have equal energy density; this is their case of 
(Jb/B = 0.6, where Ob is the dispersion of the random component and B is the strength 
of the uniform component. We compare our data with their results in the model (Fig. 
[T8h). Optical depths tk s are calculated from co lor excess E(H — Ks) using the relations 



A K JE(H - Ks) = 1.44 (INishivama et all 120061 1 and t Ks = A Ks /2.5 logi e. Some of the 
bulge sources are distributed below the boundary with <tb/B = oo. Taking errors of 5px s 
< 1% into consideration, these measurements can move to the region above the boundary. 
The bulge sources show relatively lower polarization efficiency than their best-fit result and 
almost lie bet ween an IB = 0. 6 and <J b/B = 1.2. The polarization efficiency toward the GC 
measured by iKobayashi et al.l (119831 ) as shown by the open circles is also lower than their 
best-fit result. These indicate that the energy density of the random component is higher 
than that of the uniform component of the magnetic field toward the GC. 

In the wave model, a magnetic field is described as a wave. The amplitude of the wave 
determines the e xtent to which the m agnetic field direction in each segment fluctuates along 
the line-of-sight. IJones et al.l (119921 ) fitted the model to the data, also concluding that the 
energy density of the magnetic field is in equipartition with the kinematic energy density of 
moving clouds; this is their case of Khis/Va = 1-0, where V rras is the rms motion of individual 
clouds of gas and dust attached to the magnetic field lines and Va is the Alfven speed. 
The comparison between our data and their results in the model is shown in Figure I18t >. 
The measurements below the boundary with V tids /Va = oo can also move to the region 
above the boundary if we take errors of Spx s < 1%. Compared to their best-fit result with 
K-his/Va = 1-0, the bulge sources show relatively lower polarization efficiency, most of which 
are distributed between V Tms /VA = 1-0 and V tios /Va = 1.5. This means that the turbulent 
energy density is higher than the magnetic energy density in the ISM toward the GC. 
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The comparison in either case suggests a higher magnetic field strength of the random 
component compared to that of the unifor m component between the GC an d us. Such a trend 
is also observed in the solar neighborhood (jHeileslll996t iFosalba et al.ll2002l ). The gas motions 
such as turbulence, gravitational contraction of dense clouds, expansion of H H regions, 
and supernova explosions distort the magnetic field and produce its random component. 
Th ese processes would be strongly active in the direction of the GC (see e.g., a review 
by iMorris fc Serabynl 119961 ). and contribute to the higher magnetic field strength ratio of 
the random to the uniform component toward the GC. The processes cause relatively large 
deviations of the magnetic field directions from the direction of the uniform component 
(the Galactic magnetic field). Low polarization efficiency and its spatial variation would be 
explained by superposition of the diffuse ISM and dense clouds with such deviations of the 
magnetic field directions along the line-of-sight. 



4.3. FLATTENING IN THE WAVELENGTH DEPENDENCE OF 

POLARIZATION 

In our results, ((3jh) is larger than {Phk s )\ the wave length dependence of polarization 
show s flattening from 1.2 5 to 2.14 /im. In previous studies (jWilking et al.lll980l . ll982t iNagata 
19901 ; ICreese et al.lll995l ). similar flattening is also seen. 



The value (p JH ) (2.08 ± 0.02) is close to the pow er law index a (L99 ± 0.02) of 
the wavelength dependence of extinction toward the GC ( iNishiyama et all 120061 ). However 
{/3hk s ) (1-76 ± 0.01) is clearly below a, and flattening cannot be seen in the wavelength 
dependence of extinction from the H to K$ band. At long e r wavelengths beyond 3 /im, both 



the w avelength depen dence of polarization (INagatal Il990l ; iMartin et all Il992l ; INagata et al. 



1994 ) and extinction ( Nishiyama et al. 2009bi and references therein) show flattening 



The wavelength dependence of polarization is determined by the polarizing grain prop- 
erties such as shape, size distribution, and composition. From a comparison between ob- 
servational data and theore tical models, the polarizing grain properties can be examined. 
Kim fc Martini ( 119941 . Il995al ] bj) fitted models of in finite cylindrical and spheroidal dust grains 
to modified Serkowski's law (IWhittet et al.lll992l ) and a single power law (j3 = 1.65) for A = 
1.64—5 /jm to examine the size (mass) distribution of dust grains. They obtained the most 
satisfactory result by adopting perfectly aligned oblate dust grains (axial ratio of 6:1). The 
resultant mass distribution has a peak at dust size of about 0.2 /jm and a shoulde r from 
the peak through dust size of 0.6 to 1.0 /im (see Fig. 3b of iKim & Martini Il995bl ). The 
shoulder is required to fit the infrared polarization with a power law behavior, which is the 
excess above the Serkowski's law. To explain the flattening (i.e., the excess above a power 
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law behavior), greater numbers of such large-size dust grains would be necessary. 



5. CONCLUSION 

We have made polarimetric imaging observations toward the GC in order to examine 
the efficiency and wavelength dependence of interstellar polarization at NIR. The results are 
as follows. 

1. The polarization efficiency of the ISM between the GC and us is lower than that of 
the ISM in the disk, by a factor of about three on average. 

2. The spatial variation of the polarization efficiency does not depend on the Galactic 
structure in contrast with those of color excess and degree of polarization. 

3. Position angles are almost parallel to the Galactic plane, suggesting that the magnetic 
field between the GC and us has the longitudinal configuration and connects to the toroidal 
magnetic field in the GC. 

4. The dispersions of position angles increase with decreasing the polarization efficiency. 
It is likely that the polarization efficiency is reduced by the different directions of magnetic 
fields along the line-of-sight (depolarization). 



5. The comparison of our data with the models by I Jones et all (119921 ) suggests that the 



random component has a higher strength than the uniform component of the magnetic field. 

6. The ratios of degree of polarization are Ph/pj = 0.581 ± 0.004 and Pk s /ph = 0.620 
± 0.002, which correspond to (3 JH = 2.08 ± 0.02 and HKs = 1-76 ± 0.01 for the power law 
indices of the wavelength dependence of polarization. The degree of polarization is higher 
than that expected from a single power law toward longer wavelengths (flattening from 1.25 
to 2.14 fim). The flattening i nvokes greater number s of aligned large-size dust grains in the 



mass distribution derived by lKim fc Martini (jl995bl ). 
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Table 1. Polarimetry of R CrA No. 88 



Study 


PJ £ 


M 




PH [ 


Vc] 


Oh [°] 


PK S [ 


»] 


^ s [°] 


this study 


3.83 ± 


0.17 


90 ± 1 


2.66 ± 


0.18 


91 ± 1 


1.60 ± 


0.21 


92 ± 3 


Whittet et al. (1992) 


3.87 ± 


0.06 


90 ± 1 


2.73 ± 


0.07 


92 ± 1 


1.69 ± 


0.08 a 


95 ± 1 



a Note that the observations by lWhittet et al.1 (|1992l ) were made in a non-standard K pass- 
band whose central wavelength is 2.04 fim. The value at 2.14 /im was calculated by the 
power law extrapolati on of the values measured at 1.64 and 2.04 //m in the same manner as 
Gerakines et al.l (|1995l ). 



Table 2. Wavelength Dependence of Polarization 



Ph/pj PKs/ph Pjh /3hk s 

mean 0.581 ± 0.004 0.620 ± 0.002 2.08 ± 0.02 1.76 ± 0.01 

standard deviation 0.076 0.047 0.46 0.25 

average error 0.055 0.063 0.36 0.37 
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Fig. 1. — The observe d area in this stu dy. The background is the JHKg composite image 
of the GC (Fig. 1.5 of lNishiyamall2005l ). Each square corresponds to 3 x 3 fields with a size 
of 20' x 20'. 
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Fig. 2. — Differences of p (left side) and 9 (right side) as a function of means of p in the 
J (top), H (middle), and Ks (bottom) bands based on the comparison of the same sources 
with 5p < 1% and 89 < 10° in overlapping regions. 
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Fig. 3. — J — H vs. if — i^s color-color diagram for the sources that are detected in all the 
bands and have 5p < 1% in all the bands. The thin and thick curves are the lo ci of dwarfs 
and giants, respectively. The data for 09— B9 dwarfs are from iKoornneefl (119831 ). and those 
for AO— M6 dwarfs and GO— M7 giants a re from Bessell fc Brett ( 19881 ). The arrow indicates 
a reddening vector whose slope is 1.72 ( Nishiyama et al. 20061 ). and its length corresponds 
to extinction of Ak s = 1 mag. The upper left cross denotes the average errors of colors for 
the sources. 
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Fig. 4. — Histograms of color excess E(H — Kg) for the bulge sources with 5pj < 1% (top), 
those with Spjj < 1% (middle), and those with Spx s < 1% (bottom). The means, standard 
deviations, and average errors of E(H — Kg) for the bulge sources are shown at the upper 
right of the panels. The dash-dotted lines in each panel show the means of E(H — Kg) for 
the disk sources with Sp < 1% in each band. 
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Fig. 5. — Map of E(H — Kg). Each pixel represents a mean of E(H — K$) for the bulge 
sources with 5px s < 1% in each cell with a size of 2' x 2'. 
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Fig. 6. — Ks band polarization vector map for the disk and bulge sources with 5px s < 1% 
and 50 k s < 10°. Each bar is parallel to the E- vector of polarization. The length of each bar 
is proportional to degree of polarization. Color excess for the bulge sources is also shown as 
color of the bars. The disk sources are shown by white bars. Each bar including the bar for 
10% scale of polarization degree at the upper right can be recognized by zooming the figure 
in the electronic edition. 
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Fig. 7. — Position angles 9 vs. if — Ifs colors for the disk and bulge sources with 5pj < 1% 
and 59 j < 10° (top), those with < 1% and < 10° (middle), and those with 5p^ s 
< 1% and 5(9^,, < 10° (bottom). The dashed lines in each panel represent the orientation of 
Galactic plane (~27°). The means, standard deviations, and average errors of 9 for the disk 
(left values) and bulge (right values) sources are shown at the lower right of the panels. 
The upper left crosses in each panel denote the average errors of 9 and H — K$ for the 
sources. 
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Gal. Long, offset [arcmin] from / = 0° 

Fig. 8. — Map of 6 Ks . Each pixel represents a mean of 9k s for the bulge sources with 5px s 
< 1% and 58k s < 10° in each cell with a size of 2' x 2'. The white pixels include no sources 
and means cannot be measured. 
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Fig. 9. — Degrees of polarization p vs. H — K$ colors for the disk and bulge sources with 
5pj < 1% (top), those with 5pn < 1% (middle), and those with 5px s < 1% (bottom). The 
means, standard deviations, and average errors of p for the disk (left values) and bulge 
(right values) sources are shown at the upper right of the panels. The upper left crosses in 
each panel denote the average errors of p and H — K$ for the sources. 
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Fig. 10. — Map of Pk s - Each pixel represents a mean of pk s for the bulge sources with 5pK s 
< 1% in each cell with a size of 2' x 2'. 
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Fig. 11. — Histograms of the polarization efficiency p/E(H — K$) for the bulge sources 
with 5pj < 1% (top), those with 5pn < 1% (middle), and those with 5p_R: s < 1% (bottom). 
The dashed lines in each panel are pj/E(H — K s ) = 25.0% / mag, pn/E(H — K s ) = 
14.5% / mag , and Vk ? I 'E{H — Ks ) = 9.0% / mag, which correspond to p max /E(B — V) = 
9.0% / mag ( ISerkowski et al.lll975l ). The means, standard deviations, and average errors of 
pj E(H — Ks) for the bulge sources are shown at the upper right of the panels (the means are 
also represented by arrows). The averages errors of p/E(H — K$) are derived from statistical 
and systematic errors of p and E(H — Ks). The dash-dotted lines in each panel show the 
means of p/E(H — Ks) for the disk sources with Sp < 1% in each band. 
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Fig. 12. — Map of pk s /E(H — K$). Each pixel represents a mean of pk s /E(H — K$) for 
the bulge sources with 5px s < 1% in each cell with a size of 2' x 2'. 
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Fig. 13. — Map of ct(9k 3 ). Each pixel represents a dispersion of 6k s for the bulge sources 
with 5px s < 1% and 89k s < 10° in each cell with a size of 2' x 2'. The cells including one 
or no source (s) are drawn by white pixels. 
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Fig. 14. — (pk s /E(H — Kg)} vs. <j(9k s ) for the cells including more than two sources. The 
crosses and error bars represent medians and standard deviations of (pk s I E(H — K s )) in 2° 
width bins of ct(#r- s ). 



-37- 




Fig. 15. — Spatial distribution of the sources that are detected in all the bands and have 5p 
< 1% and p > 10 5p in all the bands. The solid lines show the observed area. 
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pu m 

Fig. 16. — pj vs. Ph (top) and pu vs. {bottom) for the sources that are detected in all 
the bands and have Sp < 1% and p > 10 5p in all the bands. The upper left crosses in each 
panel denote the average errors of p. 
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Fig. 17. — Histograms of (3jh (top) and (3hk s (bottom) for the sources that are detected 
in all the bands and have Sp < 1% and p > 10 5p in all the bands. The means, standard 
deviations, and average errors of (3jh and (3hk s f° r the sources are shown in Table El 
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Fig. 18. — Degrees of polarization p Ks vs. optical depths t Ks for the bulge sources with 
8pK s < 1%- The cross indicates the means of px s and tk s for the disk sources with 5px s 
< 1%. The lines show (a) the results for five values for the para meter an IB in th e two- 
component model and (6) those for V rras /Vx in the wave model by IJones et all ( 119921 ). The 
thick lines are their best-fit results. The circles show average polarization seen toward several 
regions a t various optical depth intervals and the central regions of several normal spiral 
galaxies ( IJones et al.lll992l . and refer ences therein) , of which the open circles correspond to 
the measurements toward the GC by lKobayashi et al.l (119831 ). 



